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Abstract 

We study the associated Higgs production with Z boson at future linear colliders in the 
framework of the minimal noncommutative standard model. Using the Seiberg-Witten map, we 
calculate the production cross-section considering all orders of the noncommutative parameter 
©/ui/- We consider the effect of earth’s rotation on the orientation of with respect to the 
laboratory frame and thus on the total cross-section, it’s azimuthal distribution and rapidity 
distribution for the machine energy ranging from 0.5 TeV to 3 TeV corresponding to the non¬ 
commutative scale A > 0.5 TeV. 
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1 Introduction 

Recently, the models of extra-spatial dimension where the gravity is strong at the TeV scale nua, 
has drawn a lot of attention among the physics community. Thanks to the TeV scale new physics, 
even if some tiny amount of space-time noncommutativity is realised in nature, one can possibly 
consider the effects of this at the currently running multi-TeV energy hadron collider like LHC or 
the upcoming TeV energy electron-positron linear colliders. Much cleaner and hadronically quiet 
environment together with the complete knowledge on incoming beam momenta at the linear collider 
permits the complete reconstruction of events. Hence, one can expect a precision measurements of 
the effects from the new physics. Although the idea is rather old, lately remarkable interests arose in 
the noncommutative (NC) field theory from the pioneering work of Snyder [3]. 

The space and time coordinates at the TeV energy become operator satisfying 

= ( 1 ) 

where is an antisymmetric matrix tensor and of dimension A is the NC scale at which 

the effect of space-time noncommutativity shows up. is the anti-symmetric c-number matrix. 

There are several approaches to study the effect of space-time noncommutativity in a field theory. 
One is the Moyal-Weyl (MW) approach in which one replaces the ordinary product between two 
functions (j){x) and 'ip{x) in terms of -k (Moyal-Weyl) product defined by a formal power series expansion 
of [4l[5l[6] 

(/* 5 ){a:) = exp f{x)g{y)\y=x- ( 2 ) 

together with the ordinary integral f d^xf{x) which has the property 

J d'^x{f-kg){x) = J d'^xigi^ f){x) = j d^xf{x)g{x) (3) 

Here f{x) and g{x) are ordinary functions on and the expansion in the star product can be seen 
intuitively as an expansion of the product in its non-commutativity. Extensive works were done on 
renormalisation of the NC field theory m lain uni umig. The non-trivial phase factor (arises after 
summing over all orders of 0 in the star product defined in the Weyl-Moyal plane) which gives rise the 
UV/IR mixing problem [131 El • Detailed collider searches of space-time noncommutativity include 
the following: Hewett et al. [ISlllellTT] and others [181119] investigated e.g. e+e ^ e+e (Bhabha), 
e~e~ e~e~ (Moller), e~j ^ 77 (pair annihilation), 77 ^ and 77 ^ 77 in 

the context of NCQED. Eor a review on the NC phenomenology see EO]. 

Another approach is the Seiberg-Witten approach in which the space-time noncommutativity is 
being treated perturbatively via the Seiberg-Witten (SW) map expansion of the fields in terms of 
noncommutative parameter 0 [2ll|22l|23l|2l]. Here the gauge parameter A and the gauge field is 
expanded as 

X^{x,e) = a(a:) + e^"AW(a:; a)+ «) + •■■ (4) 

A,{x,e) = A,{x) + e>^^A^X{x) + e>^^e^^A%,^{x) + --- (5) 
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The advantage in the SW aprroach over the Weyl-Moyal approach is that it can be applied to any 
gauge theory and matter can be in an arbitrary representation. Reasonable progress has been made 
in NCQED using SW map as far as it’s quantum structure, perturbative renormalization is concerned 
[25l [26l Wf[ [28l [29l [3Ql [311 (3^ (331 |34l |35l [36] . In particularly, the UV/IR mixing arised in an arbitrary 
non-abelian noncommutative gauge theory has been studied [ISlIMlEHlEIllllllllllslESlIle]. Using 
SW expansion of the NC fields, Bichl et al. m show that the photon self-energy in NCQED 
is renormalizable to all orders in 0. Anomalies and renormalizability of SW ^-expanded NCSM are 
studied in [million]. The SW map based 6 >-exact model two-point functions for photon and neutrino 
are discussed in [sniiaiiiiiisiisei. By calculating the one loop self-energy correction to the massless 
fermion to order 0, Horvat et al. (See [33]) show that the UV/IR mixing effects can be made under 
control. In the SW map, the first several orders of the expansion can be written in a simple form 
by introducing certain generalized star products HZIEH]. Such an expansion enables one to treat all 
orders of 0 at once in each interaction vertex, thereby allows one to compute nonperturbative results. 
In [31] , this technique was used to compute the fermion one loop correction to the photon two point 
function of a NCQED model using SW fields. 

Using this SW map Calmet et al. [39] |40] first constructed the minimal version of the noncom¬ 
mutative standard model (mNCSM in brief). They derived the 0(0) Eeynman rules of the stan¬ 
dard model interactions and found several new interactions which are not present in the standard 
model. All the above analyses were limited to the leading order in 0. Das et al. first analysed the 
e+e“ ^ 7 , Z ^ to order 0^ (without considering the effect of earth’s rotation) [41] . There 

exists another version: the non-minimal version of the NCSM (nmNCSM in brief) where the triple 
neutral gauge boson coupling arises (absent in the mNCSM) naturally in the gauge sector. This 
model (nmNCSM) was first formulated by Melic et al. [42j|43]. Interesting phenomenological studies 
comprising triple gauge interaction are available in the literature [HIlSlIlllllTlIlH]. 

The direct test of space-time noncommutativity comes from the two decays Z ^ 77 , gg (forbidden 
in the SM at the tree level). Using the experimental bound < 1-3 x 10“^ GeV and < 

1.0 X 10“^ GeV, Behr et al. [44] shows the bound on the NG scale ^ 1 TeV. Taking the SM fields 
in the enveloping algebra, Calmet et al. gaisn] shows that the bound on the NC scale ^10 TeV, a 
rather weak bound. 

The impact of neutrino-photon interaction in the noncommutative space-time on the cooling of 
stars [51], on the primordial nucleosynthesis and ultra-high energy cosmic ray [52] [53] have been 
studied in detail. Assuming the plasmon decay to a pair of neutrino may contribute substantially to 
the star energy loss, the author P. Schupp et al. EB obtain a lower bound A > 81 GeV. The non¬ 
observation of large neutrino-nucleon cross-section for ultra high energy neutrinos (of energy 10 ^^ GeV 
at neutrino observatories gives rise a lower bound on A which is as high as 900 TeV [53]. The invisible 
Z boson decay in covariant theta-exact NGSM was studied in [54], The impact of NG space-time on 
Quarkonia decay into two photons [55l [56] and K ^ m decay (forbidden in the SM) have been 
investigated in detail. However, the experimental upper bound on such rare decays are too weak to 
obtain any bound on A. Early collider searches of space-time noncommutativity include the work by 
Kamoshita et al. [58] [59] . Das et al. (one of the current authors) studied in detail the Bhabha and 
the Moller scattering [ 6 Q[[ 6 T| . muon pair production [62][63] in the non-minimal NGSM scenario. The 
non-commutative parameter can be fundamental constant in nature that has a fixed direction 

in the celestial sphere. Hence, the daily modulation effect of earth’s rotation can be observed in the 
noncommutative phenomenology. Now after the finding of ~ 125 GeV scalar boson at the Large 
Hadron Collider, CERN in 2012 , 7 TeV and 8 TeV LHC data has also established that this scalar is 
most likely be SM Higgs which is the last peace of SM to be discovered. Now the most important 
questions are as follows: whether there is any new physics beyond the standard model at the multi 
TeV scale which can show up in present and future colliders? And, whether the Higgs boson can throw 
some light in this? At the TeV energy scale where the space-time becomes noncommutative(NC), 
whether such a NC space-time can play an important role in Higgs boson production and it’s decay? 
In an earlier work one of the current authors studied the Higgs boson pair production in the nmNCSM 
scenario [64] [65] . More recently, Wang et al. [ 66 ] studied the Higgs-strahlung process in the context of 
noncommutative standard model at linear collider without considering the effect of earth’s rotation 
into account. Besides the Moyal-Weyl product and Seiberg-Witten map, there is another way to 
deal with space-time noncommutativity. In this approach, one don’t use any deformed products like 
the Moyal product: one map the field theory model onto a kind of matrix model by representing the 
associative algebra modeling the noncommutative space as an algebra of operator acting on a suitable 
Hilbert space. 
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In this paper, we explore the associated Higgs (of mass 125 GeV) production with Z boson at 
the TeV energy Linear Collider in the context of non-minimal standard model. We studied the 
time-averaged cross-section, azimuthal distribution, NC correction to the cross-section, the rapidity 
distribution (of the final state Z and H boson) in detail considering the effect of earth rotation into 
account and explore the important role played by the orientation of the NC vector (0) on different 
observables which can be tested in the upcoming TeV energy linear collider. We organise the content 
as follows. In Sec. we discuss the noncommutative vertices involved in the process and obtain the 

squared matrix amplitude for the Higgs-strahlung process e“e+ —> ZH and obtained the expression 
of cross-section and angular distribution with and without considering the effect of earth rotation 
into our analysis. 

Sec. is devoted for the numerical analysis. We investigate the impact of space-time noncommu¬ 
tativity on cross-section, angular distribution, rapidity distribution. We discuss about the potential 
relevance of the TeV scale noncommutative geometry at the linear collider. Finally, we summarise 
our results and conclude in Sec. O 

2 e“e+ ^ ZH process in nmNCSM 

Associated Higgs production with Z boson and the vector boson fusion through W oi Z bosons are 
two most significant Higgs production channels in linear collider. Both of these processes involve the 
VVH {y = Z) couplings which relates the SM symmetry breaking and precise measurements 

and firm constrains on anomalous couplings can be expected from linear colliders. For our present 
calculation, the tree level Feynman diagram for the Higgs-strahlung process which is a s-channel 
process is shown below in Fig.[2 We use the notation for the incoming electron and positron momenta 



Z^ 

Figure 1: Representative Feynman diagram for the proeess e“(pi)e+(p 2 ) ^ Z{ps)H{p 4 ) both in SM 
and NCSM with different strueture of eouplings. As ean he noted from the vertiees given in the text, 
in the very large NCSM seale (A ^ oc) extra tensor struetures disappear to reproduee SM eouplings. 


as Pi and p 2 , whereas the outgoing Z boson and Higgs(iL) boson momenta are ps and p 4 , respectively. 
In terms of the noncommutative parameter we use the Feynman rule [66] for e~ —e'^ — Z vertex 
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and for Z — Z — H vertex 
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Neglecting the masses of incoming electrons and positrons, we consider the vertices above contain all 
orders of 0 (oc terms. Using these Feynman rules, we find the squared amplitude (spin-averaged) 
as, 

where the quantity psQp^ (the argument of cosine function appeared in above) is given by 


(8) 


(^) as Vs » Mz,Mh (9) 

where f{0y) = cosO + sinO {sincf F coscj)). Here 6>, (f are the polar and azimuthal angles of the 
outgoing Z boson. A ( 5 , (the Kallen function), in terms of Higgs mass (Mh) and Z boson 
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mass {Mz) are given by 


A(s, M|, M^) = s^ 2sMl - 2sMfj - 2M|M^ ^ 5 ^ as > Mz, Mr 


It is important to note that the oscillatory behaviour that the cross section and other angular depen¬ 
dence of the azimuthal distribution (as we will see later) is due to the presence of the function /(^, 0) 
defined above. 

Note that the SM squared amplitude term gets recovered in limA^oo \^\‘^ ncsm equal to 
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( 10 ) 


Since the noncommutative parameter is considered as fundamental constant in nature, it’s di¬ 
rection is fixed with respect to an inertial (non rotating) coordinate system. Now the experiment is 
done in the laboratory coordinate system which is located on the surface of the earth and is moving 
by the earth’s rotation. As a result 0^^^ (©£;, ©b), fixed in the primary co-ordinate system, will also 
vary with time in the laboratory frame and this must be taken into account while making any serious 
phenomenological investigations of space-time. 

To study the effect of earth’s rotation in the non-commutative space-time, we follow the notation 
by Kamoshita [58]. Let ix^ Jy and kz be the orthonormal basis of the primary (non rotating) 



Figure 2: In the figure the primary eoordinate system (X-Y-Z) with the Z-axis along the earth axis of 
rotation and the laboratory eoordinate system (i — j — k) for an experiment on the earth are shown. 
The generie NC veetor 0 is shown in the X- Y-Z system where g and f eorresponds to the polar and 
the azimuthal angle. In the above time varying ( = ujt where uj is a eonstant. (6, a) defines the 
loeation of the laboratory. 

coordinate system (X-Y-Z). In the laboratory coordinate system {x — y — z), the bases vectors are i, j 
and k. The base vectors of the primary (non rotating) coordinate system can be written in terms of 
the laboratory coordinate system as [58] 

ix (Ca<S(^ T <S(5 Sq,C(^) i T C(5C(^ j “h (^S(iS(^ S^C(iC(^^ k^ 

jx = (-CaC^ + SsSaS(^) i + j + (-5^0^ - SsCaS^) k, 

kx = -CsSa i Y Ss j -I CsCa k (11) 

Here Ca = cosa, Sa = sma, cs = cosS, C(^ = cos( etc. In Fig. [^we have shown the primary(X —Y —Z) 
and laboratory(i—jf — /c) coordinate system. Note that the primary Z axis lies along the axis of earth’s 
rotation and {6, a) defines the location and orientation of e~ — e+ experiment on the earth, with 
— 7 r /2 < S < IT/2 and 0 < a < 27r. Because of earth’s rotation the angle ( (see Fig. 2 ) increases with 
time and the detector comes to its original position after a cycle of one complete day, one can define 
C = ojt with angular velocity uj = 27rjT^ay and T^ay = 23h56m4.090535. In Fig. the orientation of 
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the NC vector 0 is described by the angles (77, C) l^he primary coordinate system with 0 < 77 < tt 
and 0 < ^ < 27r. The differential cross-section of the e+e“ ZH scattering is given by 


da 

d cos Odcj) 


647r^5^ 


NCSM 


( 12 ) 


and the total cross-section 


PIT n27T 7 

a= de d4>- - - (13) 

Jo Jo d cos 9 d(p 

In above, the spin-averaged squared-amplitude I^Patcsm given by Eqnj^ To extract the effect of 
this lab-rotation coming through noncommutative effect, one takes the average of the cross-section a 
or it’s distributions over the sidereal day Tday Here, we consider (cr)^ = ^ 

[SHI [63l [65] . Here a = cr(y^. A, 77, t). The cross-section is calculated using the center of mass frame 

of the e+e” —5- ZH process in which 4 momenta of the incoming and outgoing particles are given 
by: 


Pi = Pe- =^(1,0,0,1), P2=Pe+ = ^( 1,0,0,-1), 

P3 = Pz = { Ez, k'sinO cos0, k'sinO sin0, k'cosO), 

Pa = Ph = { Eh, — k'sinO cos0, — k'sinO sin0, —k'cosO). 


where k' = ^ y 7 and 0 is the scattering angle made by the 3-momentum vector ps of Z boson with 

the k axis (the 3-momentum direction of the incoming electron e~) and (p is the azimuthal angle. The 
time dependence in the cross-section or it’s distribution enters through the NC parameter 0(= Qe) 
which changes with the change in ( = ujt. The angle parameter ^ appears in the expression of 0 
through cos{ujt — or sin(cc;t — <^) as the initial phase for time evolution gets disappeared in the time 
averaged observables. So one can deduce 0 i.e. the NC scale A and the orientation angle 77 from the 
time-averaged observables. 


3 Noncommutative effects on Cross-section and angular dis¬ 
tributions 

We analyse the Higgs-strahlung process, e“(pi) e+(p2) ^ E{ps) H{p 4 ) in presence of the non¬ 
commutative standard model in TeV energy linear colliders. We assume both the final Z boson and 
Higgs boson are produced on-shell and reconstructed from their decay products which is not prob¬ 
lematic in a linear collider by choosing suitable decay channels. The initial unpolarised electron and 
positron beams are colliding back to back with half of machine energy without considering any effects 
from ISR whereas the final particle momenta can be defined in terms of polar (0) and azimuthal 
angle (0) as they are produced inside the collider. We further study this process with and without 
taking into consideration the effect of earth rotation into noncommutative space-time. In this context, 
we probe the non-commutative scale A for the machine energy ranging from 500 GeV to 3000 GeV 
from the associated production of Higgs with Z boson. Here we utilise the total cross-section rate, 
azimuthal distribution and rapidity distributions and their time-averaged estimates to discriminate 
the new physics. In our analysis we have used the Higgs boson mass as 125 GeV and we set the 
laboratory coordinate system by taking {6, a) = (7r/4,7r/4) which is the OPAL experiment at LEP. 
In Eig. [^we display the total cross-section as a function of the center of mass energy energy ^/s for 
different values of non-commutative scale A. The uppermost curve in this plot refers the expected 
contribution from SM production (which corresponds to the NGSM value at the limit of NC scale 
A ^ 00 as we already noted down in previous section). While going below in the same plot, the 
curve next to the topmost one corresponds to the scale A = 1.5 TeV and so on. On the other hand, 
the lowermost curve corresponds to A = 0.6 TeV. The deviation from the SM plot starts getting 
manifested at and above = 1 TeV. Eor example, at the machine energy =1.5 TeV, we see 
that the overall cross-section gradually increases with the increase in the scale A from 0.6 TeV to 
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Figure 3: The total cross-section a (in fb) for associated production of Higgs with Z boson is shown 
as a function of the machine energy y/s (in GcV). Different lines represent the choice of different 
non-commutative scale A which is ranging from 0.6 TeV to 1.5 TeV. The topmost curve (solid line) 
corresponds to the expected Standard Model production cross-section which is essentially NC cross- 
section at the limit A ^ oo. 



yfs (in GeV) 



Figure 4: On the left, (a) noncommutative correction to the total cross-section Act (in fb) for asso¬ 
ciated production of Higgs with Z boson is shown as a function of the linear collider machine energy 
^/s (in GeV) for different A ranging from 0.6 TeV to 1.5 TeV. On the right, (b) the same quantity 
is plotted as a function of NO scale A for different machine energy. 


1.5 TeV before merging with the SM value for larger NC scale. Expectedly, maximum deviation from 
that of SM in observed for the lower values of A. Also, note that for a given machine energy the NC 
cross-section is always less than that the corresponding SM value. That is simply followed from the 
Eq-i Moreover, one can notice that at higher machine energy the NC cross-sections are not always 
simply falling with monotonous regularity. The ripple effect appears and become prominent at some 
higher machine energy subject to each NC curves. This is an effect of tensor structure coming into 
the Eq. In the next subsection we would explore this as a characteristic feature from NC effects. 


3.1 Noncommutative correction 

To explore and quantify the noncommutative effect manifested in production cross-section and orig¬ 
inated from the tensorial strunture as in we define the NC correction with respect to the SM 
value Act as. 


N(7 = (Tnc - CFSM’ (14) 

In Fig. I^a), we have demonstrated the variation of Act as a function of the machine energy y/s 
for different values of A. Now, we can spot it very clearly that Act first decreases (becomes more 
negative as cfnc < (^sm) with the increase in y/s and reaches a minimum for its first trough 0 Note 
that the first minimum occurs at ^/s = 2 x (3)^/^ x A (See Eqnj^= 1579(2105) GeV corresponding 
to A = 600(800) GeV,whereas the second minimum occurs at ^/s ^ 2700(3600) GeV for the same 
A. After that they decreases with the increase in ^/s and pass through an oscillatory phase and 

^However additional subsequent troughs remain subdominant to be observed in this figure and would be clearly 
prominent in our next figure where we would consider the relative correction 
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Table 1: The NC correction Sr against the NC scale A (in GeV) is shown corresponding to different 
machine energy ^/s. Primary peak values and corresponding NC scales are shown in bold. 


ys = 500 GeV 

= 1 TeV 

Vs = 1.5 TeV 

Vs = 3 TeV 

A (GeV) 

Sr 

A (GeV) 

Sr 

A (GeV) 

Sr 

A (GeV) 

Sr 

100 

-0.543 

200 

-0.515 

300 

-0.508 

600 

-0.504 

158 

-0.608 

300 

-0.569 

450 

-0.564 

700 

-0.485 

160 

-0.607 

330 

-0.608 

497 

-0.608 

800 

-0.452 

250 

-0.229 

400 

-0.485 

600 

-0.492 

999 

-0.608 

300 

-0.120 

600 

-0.139 



1100 

-0.572 

600 

-0.008 








eventually becomes a ffat curve asymptotically meeting at the line Aa = 0, which is the case for large 
A limit (the SM limit). As for example, for A = 0.6 TeV, we find Aa = —2.7 fb at ^/s = 1.5 TeV 
and for A = 0.6 TeV, Aa = —0.4 fb (minimum) at ^/s = 4.0 TeV. As we vary A = 0.6 TeV to 
A = 1.5 TeV, we see that the height/depth of the trough Aamin decreases and it’s location (the 
value of ^/s) also gets changed. In Fig. |^b), we have plotted the variation of the difference Act as a 
function of the NC scale A corresponding to the different machine energy ^/s = 0.5, 1.0, 1.2, 1.5 and 
3 TeV, respectively. We see that the deviation decreases with the increase in A for a fixed machine 
energy and for large A the NC effect gets disappeared as expected. For ^/s = 1 TeV, |Acr| changes 
by an amount 6 fb as A changes from 400 GeV to 1500 GeV. 

As noted earlier, we next define the relative correction of this NC cross-section by normalising the 
noncommutative correction to the total cross-section by the SM value: 

5r = —. (15) 

O'SM 


In Fig. we have shown the variation of Sr against the machine energy ^/s for different values 
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Figure 5: On the left, (a) the ratio Sr is plotted as a funetion of ^/s (in GeV) for different A values. 
On the right, (b) Sr is shown as a funetion of A (in TeV) for different maehine energy. 

of A (left figure), whereas on the right figure, we have shown the dependence of Sr against the 
NC scale A for few different machine energy ^/s. From the left figure, we see that Sr becomes 
maximum at its first peak when ^/s = 1.5, 1.8, 2.1, 2.4 ,2.7 TeV and 3.0 TeV corresponding to 
A = 0.5, 0.6, 0.7, 0.8, 0.9 TeV and 1.0 TeV. For A = 0.5 TeV at = 1.5 TeV, the deviation Sr 
reaches it’s minimum value —0.6085. This is called primary minimum (crest). Similarly the secondary 
minimum (crest) corresponds to Sr = —0.5456 and the next to it corresponds to Sr = —0.435. On 
the right plot, we have shown how Sr varies with A as the machine energy increases from 0.5 TeV to 
3 TeV. For A < 1 TeV, we see that for a wide range of machine energies, Sr converges to —0.5. In 
Table we have displayed the value Sr for different A corresponding to different machine energies. 
We also listed several values of Sr for different A corresponding to CLIC energy (y^ = 3.0 TeV). We 
find that Sr(= —0.608596) is maximum corresponding to A ~ 1 TeV. 
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3.2 Angular distributions in absence of earth rotation 

The angular distribution of the final state scattered particles is a useful tool to understand the nature 
of new physics. Since the noncommutativity of space-time breaks Lorentz invariance including rota¬ 
tional invariance around the beam axis, this will lead to an anisotropy in the azimuthal distribution 
of the cross-section i.e. the distribution will depends strongly on (j). In the standard model the az¬ 
imuthal distribution for one of the final particles is found to be flat. However, in the noncommutative 
standard model(NCSM) due to the presence of the tensor 6>-weighted dot product i.e. terms like 

PsOpA (cosOsinO {sincj)coscj))) (Eqn. qL these distributions are no longer remain flat. 

We see (Eqs.[^[^ that the squared amplitude is a oscillatory function of (6>, (/>) and is further distorted 


dtr/d^ (fb/radian) 



Eigure 6: The azimuthal distribition ^ (in fb/rad) is plotted as a funetion of (j) (in rad) for different 
A values. Displayed are four plots (a), (b) and (e) eorresponding to different maehine energy ^/s = 
0.5, 1.0, 1.5 and 3.0 TeV are shown. 

by an oscillatory function of (y^. A) (of the form ^ cos | {cosO + sinO {sincf + co50))| ). In 

Fig.§ we have shown the azimuthal (j) distributions for different A values corresponding to the machine 
energy ^/s = 0.5, 1.0, 1.5 and 3.0 TeV, respectively. The distribution has several maxima(minima) 
located at 0 = 2.4, 5.4 rad (0.8, 4.0 rad). Eor a given machine energy (e.g. ^/s = 1.0 TeV), the 
height of the peaks decreases with the increase in A. We see that as the machine energy is increased 
to 3 TeV, the peaks corresponding to A = 0.5, 0.6, 0.7 TeV located at 0 = 2.4, 5.4 rad gets smeared, 
while the peaks corresponding to A = 1.0, 1.5 TeV still survives. The fluctuation observed at the 
crest and trough of the figure corresponding to machine energy ^/s = 3 TeV is due to the fact that 
^ is a oscillatory function of (j) and as well as the oscillatory function of (y^. A). At higher energy 
^ is dominated by y^ and A. 


Next, we define the rapidity of a particle as 


y=hn(^ —^ 
^ 2 \E + P^ 


where E, the energy and P^, the 2;-component momentum of the particle {Z boson or the Higgs boson 
H). In Eig.[^ we have plotted the distribution dafdy as a function of rapidity y for different cases of 
machine energies. In the leftmost Eigure where the machine energy is fixed at y^ = 0.5 TeV, we see 
that the height of the peak (located at ^ = 0) increases with the increase in A. The topmost curve 
corresponds to the SM curve (A ^ oo). As we move towards left (i.e. increase the machine energy 
from 0.5 TeV to 3.0 TeV, the height of a peak (corresponding to a particular A value ) decreases 
and gets flattened. When the machine energy is equal to 3 TeV, the peak at ^ = 0 corresponding to 
A = 0.6 TeV split into two. Note that the rapidity distribution of Z boson or Higgs particle are the 
same. 
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Figure 7: The rapidity distribution ^ (in fb) is plotted as a funetion of the rapidity y for different A 
values. Three plots (a), (b) and (e) are shown eorresponding to maehine energy y/s = 0.5, 1.0 and 
3.0 TeV; respeetively. 


4 Consequence of earth rotation on the cross-section and an¬ 
gular distributions 


4.1 Cross-section, diurnal motion in presence of earth rotation 


It was pointed out earlier that one significant aspect for noncommutative effect can be originated 
from directionality of fundamental NC fix points. The experiment is done in the laboratory frame 
attached to the earth surface which is rotating, whereas directions related to the NC parameter 0^^^ 
has a fixed direction in the celestial sphere. In Sec.2 we described the notation to parametrise the 
processes in rotating frame. This rotation can have a direct but subleading impact in daily modulation 
on the inherent structure of the interaction couplings. Noncommutative contributions and angular 
dependence due to additional tensorial were presented in Se. 3.1, where this effect was not considered. 
Now we would like to analyse the effect of earth rotation on the orientation of the NC vector 0^; and 
thus on the cross-section and angular distribution of the associated Higgs production. Since the cross- 
section and the angular distributions are function of time, we made a time-averaged ( i.e. averaged 
over the side-real day T^ay) estimate of the total cross-section, correction, azimuthal distribution 
and rapidity distributions to account for this additional effect coming from the new physics. The 
laboratory coordinate system is being set at {6, a) = (7r/4,7r/4) which is the OPAL experiment at 
LEP. Our choice of same lab-system enables one to directly compare the consequence due to the earth 
rotation. The time-averaged azimuthal distribution and the cross-section are defined as. 


da \ 

d(f> / rj. 


1 

Tday Jo d(f) Tday 


I rTday /*! 

day Jo J —1 


da 


_i dcosO d(f 
da 


dcosO dt. 


1 rTday 1 pTday /*27r 7 

Mt = / (^dt = / / / -- -dcosO dcf dt 

Tday Jo Tday Jo J-1 Jo d<f dcOsO 


(16) 

(17) 


where Tday = 23h56m4.09053s, the sidereal day. 

In Fig. we have shown the time-averaged cross-section {a)T as a function of the machine energy 
^/s corresponding to 77 = 0, 7r/4 and 7r/2. Only one NC scale (A = 0.6 TeV) is chosen for demon¬ 
stration comparing with the corresponding non-rotation plot in Fig. 3. The plots corresponding to 
77 = 0, TT and 77 = 7r/2 are seen to be nearly overlapping with a narrow effect due to different choices 
of 77 values. To demonstrate the variation due to this parameter, we next define the NC-correction 
to the cross-section (time-averaged) as {Aa)T = {cfnc)t — {(7sm)t‘ In Fig. 8(b), we have plotted 
{Aa)T as a function of 77 corresponding to A = 0.6, 0.7, 0.8, 0.9, 1.0 and 1.5 TeV for a fixed machine 
energy ^/s = 1.5 TeV. The plot shows a peak at 77 = 7r/2 which corresponds to the fact that { aNc)T 
is larger at that value irrespective to the A chosen. However this maximum deviation is quite small; 
around 0.1 fb corresponding to A = 0.6 TeV and even less for larger A. 

In Figj^ and Figj^, we have shown {Aa)T and {Sr)T as a function of ^/s for A = 0.6 TeV and 
77 = 0, 7r/4 and 7r/2, respectively. We see that the curves corresponding to different 77 are almost 
overlapping for different machine energy values except near the region ^/s = 1.5 TeV. From the 
lowermost to the topmost curves 77 corresponds to 0, 7r/4 and 7r/2, respectively. We see that for 
77 = 7r/2, the deviation (Acr)^ and the normalised deviation = {Aa)T/{ cfsm)t are minimum, 

yielding { aNc)T is largest for 77 = 7r/2. We set 77 = 7r/2 in the rest of our analysis. As expected from 
these results, one finds that the variation of time-averaged quantity {Aa)T and {6r)T with respect to 
machine energy ^/s (for different A) or A (for different ^/s) remains very similar to the plots shown 
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- A=0.7TeV 


. A=0.8TeV 

- A-0.9TeV 



- A=1.0TeV 
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Figure 8 : (a) The time-averaged total eross-seetion {cf)t fb) for assoeiated produetion of Higgs with 
Z boson is shown as a funetion of the linear eollider renter of mass energy ^/s. The topmost eurve 
(solid line) eorresponds to the expeeted Standard Model produetion eross-seetion whieh is essentially 
NC eross-seetion at the limit A ^ oo. The three other plots (below the SM plot) eorresponds to 
77 = 0, IT/A and it( 2 and A = 0.6 TeV are found to be almost overlapping, (b) The time-averaged 
NC eorreetion to the eross-seetion (Acr)^ (in fb) is shown as a funetion of orientation angle rj of 
the NC veetor for a fixed maehine energy y/s = 1.5 TeV. The different plots eorrespond to A = 
0.6, 0.7, 0.8, 0.9, 1.0 and 1.5 TeV. 



Figure 9: In (a) the time-averaged NC eorreetion to the eross-seetion (Acr)^ (in fb) is shown as a 
funetion of the maehine energy ^/s (in GeV) for A = 0.6 TeV. The different plots eorrespond to 
77 = 0 , 7 r/ 4 , 7 r/ 2 . In (b), the normalised eorreetion {Sr)T 'is shown as a funetion of y/s for the same 
above set off] values. 


in Figand Fig. in absence of the earth rotation with a subleasing shift. We skip the repetition of 
these plots for the brevity. Magnitude of the shifts for different parameters are better expressed in 
a representative table similar to the one we discussed in last section. In Table we have displayed 
time-averaged {Sr)T for different A corresponding to different machine energies and r] = irf2 after 
taking consideration of earth rotation effect. Note, for example, the shift in the ratio |((^r)T| = 0.604 
form our non-rotating estimate of \{Sr) \ = 0.608 for values of A = 158 GeV and ^/s = 500 GeV. 

Diurnal modulation in the production signal in our lab frame can also appear as a distinctive 
feature of fundamental NG fixed points. This modulation would depend upon 77 together with the 
NC scale A and machine energy ^/s, although the phase of these oscillation are fixed by the choice of 
experiment location. 

To take into account this modulation, we define the quantity A(t) = a{t) — {a{t))T- This is the 
deviation of the Higgs-strahlung production cross-section a{t) at any time from the time-averaged 
cross-section {a{t))T over the period of sidereal day, Tday = 23.934 hours. In Fig. we have plotted 
A(t) as a function time fraction for the machine energy ^/s = 0.5, 1.0, 1.5 and 3.0 TeV, respectively. 
In each plot, the NC scale is chosen as A = 0.6, 0.8 and 1.0 TeV and the orientation angle 77 = 7t/2. 
From the figures corresponding to ^/s = 0.5, 1.0, 1.5 TeV, we see that at a given machine energy, as 
we increase A from 0.6 TeV to 1.0 TeV, the fluctuation gets diminished and eventually it becomes 
zero in the A ^ 00 limit (the SM result). The plots show peaks at t = 0.35 Tday and 0.85 Tday 
times of the day, where Tday = 23.934 hours. For the machine energy ^/s = 0.5, 1.0, 1.5, we see that 
with the increase although the location of several peaks/dips remains unchanged and the fluctuation 
pattern remains almost same, however it’s magnitude at any particular point of time in a day for 
a give A changes largely. However, for ^/s = 3.0 TeV, we see something different behaviour: there 
are dips and peaks in the plot corresponding to A = 0.6 TeV and 0.8 TeV at some Tratio^ the same 
Tratio corresponds to dips and peaks for the A = 1.0 TeV plot i.e. they are out-of-phase. Also 
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Table 2 : The NC correction {Sr)T against the NC scale A is shown corresponding to different machine 
energy ^/s and orientation angle of the NC vector 77 = 7 r/ 2 . Primary peak value and the corresponding 
NC scales are shown in bold. 


^ = 500GeV 

^ = lOOOGeV 

Vs = 1500GeV 

Vs = SOOOGeV 

A GeV 

{6r)T 


A GeV 

{6r)T 

A GeV 

{6r)T 

A GeV 

{6r)T 

600 

700 

800 

999 

1100 

-0.504 

-0.485 

-0.456 

-0.602 

-0.563 

100 

158 

160 

250 

300 

600 

-0.541 

-0.604 

-0.603 

-0.226 

-0.118 

-0.008 

300 

450 

497 

600 

-0.508 

-0.563 

-0.602 

-0.483 

200 

300 

330 

400 

600 

-0.514 

-0.568 

-0.602 

-0.476 

-0.136 



(c) 


rratio = t(in hrs)/rday(in hrs) 


Tratio = t(in hrs)/rday(in hrs) 


Figure 10: The diurnal modulation A{t) in the produetion signal is plotted as a funetion of time 
fraetion of sidereal day Tratio{= t/Tday) for the maehine energy y/s = 0.5, 1.0, 1.5 and 3.0 TeV, 
respeetively. The NC seale is ehosen as A = 0.6, 0.8 and 1.0 TeV and 77 = 7 r/ 2 . 


interestingly, the height of the peak/dip decreases with the decrease in A, contrary to the one found 
in y/s = 0.5, 1.0 and 1.5 TeV cases. 

4.2 Angular distributions in presence of earth rotation 

The anisotropy emarged from the breaking of Lorentz invariance as mentioned earlier, persists in the 
time-averaged (averaged over the side-real day T^) azimuthal distribution of the cross-section (^)t 
can act as a signature of space-time noncommutativity which is found to be absent in many theories 
beyond the Standard Model physics. 

In Fig. we have shown the time averaged azimuthal distribution of the cross-section (^)t 
corresponding to different values of NC scale A = 0.5, 0.6, 0.7, 1.0 and 1.5 TeV. The machine 
energy ^/s is fixed at 0.5, 1.0, 1.5 and 3.0 TeV, respectively. Each distribution has maxima(crest) 
and minima(trough) located at 0 = 2 . 2 , 5.4 rad ( 0 . 8 , 3.8 rad). At a fixed machine energy (say 
^/s = 1.0 TeV), the height of the peaks decreases with the increase in A. For 77 = 7 r/ 2 , as the machine 
energy is increased from 0.5 TeV to 3 TeV, the peaks corresponding toA = 0.5, 0.6, 0.7 TeV located 
at 0 = 2 . 2 , 5.4 rad gets smeared, while the peaks corresponding to A = 1.0, 1.5 TeV still survives. 
The fiuctuation observed at the crest and trough of the figure corresponding to machine energy 
^/s = 3 TeV is due to the fact that (^)t is a oscillatory function of (j) and as well as the oscillatory 
function of (y^. A). As we see that the behaviour of (^)t at large energy {\fs) is dominated by the 
machine energy at a given A. 
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Figure 11: The time-averaged azimuthal distribution of the eross-seetion (^)t (fh/rad) is shown as a 
funetion of the azimuthal angle (f (in radian) forrj = 7 r /2 and the maehine energy y/s = 0.5, 1.0, 1.5 
and 3.0 TeV (whieh eorresponds to hqures a,b,e and d). The different plot in eaeh figure eorrespond 
toA = 0.b, 0.6, 0.7, 1.0 and 1.5 TeV. 


In Fig 


■13 

llstc 


we have also shown the time averaged rapidity distribution (^)t against the rapidity y 


of the finalstate particle for different value of the machine energy and p = 7 r/ 2 . On comparing Fig. 12 


dtr/dy(in fb) do-/dy(in fb) d<r/dy(in fb) 



Figure 12: The rapidity distribution (^)t (i'^ GeV) is shown plotted as a funetion of the rapidity y 
for T] = nf2 and the maehine energy ^/s = 0.5, 1.0 and 3.0 TeV (whieh eorresponds to figures a,b and 
e). The different plot in eaeh figure eorrespond to A = 0.5, 0.6, 0.7, 1.0 and 1.5 TeV. 

and Fig. we see that the behaviour of the rapidity distributions in both cases (in the presence and 
absence of the effect due to earth rotation) are almost identical, although the magnitude differs from 
each other slightly. 

5 Summary and Conclusion 

In this paper, we have investigated the associated Higgs production with Z boson at the future TeV 
energy linear collider. We did our calculation in the framework of the non-minimal noncommutative 
standard model(nmNCSM) using the Feynman rules involving all orders of the noncommutative 
parameters 0^^^ with (or without) considering the effect of earth rotation. We found that the total 
cross-section cr(e“e+ ^ ZH) departs significantly from the standard model value as the machine 
energy starts getting larger than 1.0 TeV with the NC contribution found to be lower than the SM one. 
Considering the effect of earth rotation in our analysis, we find from (Acr)^ vs rj plot that for 77 = 7 r/ 2 , 
(Acr)^ becomes minimum (yielding {(Jnc)t to maximum for A ^ 0.6 TeV. The time-averaged NC 
correction {/A(t)t and the relative correction {5r)T is found to be minimum (a trough) (yielding 
{(^nc)t to maximum) at ^/s = 1.5 TeV for the NC scale A = 0.5 TeV. The trough depth decreases 
with the increase in A and becomes zero as A ^ 00 (the SM value). The diurnal modulation of the NC 
signal is found to be quite interesting. We plot cr(t) —(cr)T is plotted as a function of tjT^ay and is found 
to have an oscillatory behaviour. At a given energy, the amplitude of oscillation gets damped with 
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the increase in A and finally it becomes zero in the limit A ^ oc (the SM limit). The time-averaged 
azimuthal distribution (^)t against (j) for different A at different machine energies is found to have 

a oscillatory behaviour because of the additional terms psOp^ ^ (^) + sinO {sincj) + coscj))) 

which content the tensorial effect from noncommutativity. Note that the distribution is completely 
flat in the standard model. The distribution shows peak at certain 0 values (similar peaks are observed 
in the case of no earth rotation) corresponding to different A at a machine energy {y/s = 0.5 TeV to 
3.0 TeV), which can be looked at in the linear collider experiment and thus test the idea of space-time 
noncommutativity in near future. 

Acknowledgments 

The work of P.K.Das is supported in parts by the CSIR Project (Ref. No. 03(1244) 12/EMR-II) and 
the BRNS Project (Ref. No.2011/37P/08/BRNS). The authors would like to thank Mr. Atanu Guha 
for useful discussions. 

References 

[1] N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B429, 263 (1998). 

[2] I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B436, 257 (1998). 

[3] H.S. Snyder, Phys. Rev. D 71, 38 (1947). 

[4] 1. F. Riad and M. M. Sheikh-Jabbari, J. High Energy Physics 08, 45 (2000). 

[5] B. Jurco and P. Schupp, Eur. Phys. J. G 14, 367 (2000). 

[6] J. M. Gracia-Bondia and J. C.Varilly, J.Math.Phys. 35, 3340 (1994). 

[7] T. Filk, Phys. Lett. B376, 53 (1996). 

[8] J. C. Varilly and J. M. Gracia-Bondia, Int. J. Mod. Phys. A14, 1305 (1999). 

[9] M. Chaichian, A. Demichev and P. Presnajder, Nucl. Phys. B567, 360 (2000),hep-th/9812180, 

[10] M. Chaichian, A. Demichev and P. Presnajder, J. Math. Phys. 41, 1647 (2000),hep-th/9904132, 

[11] M. Van Raamsdonk and N. Seiberg, J. High Energy Physics 03, 035 (2000). 

[12] 1. Ya. Aref’eva, D. M. Belov, and A. S. Koshelev, Phys. Lett. B476, 31 (2000). 

[13] S. Minwalla, M. Van. Raamsdonk and N. Seiberg, J. High Energy Physics 02, 020 (2000). 

[14] M. Van Raamsdonk, J. High Energy Physics 11, 006 (2001). 

[15] J. Hewett, F. J. Petriello and T. G. Rizzo, in Signals for Noncommutative QEC at High Energy 
e+e- colliders, eConf C010630, E3064 (2001). 

[16] J. Hewett et al. , Phys. Rev. D 64, 075012 (2001). 

[17] J. Hewett et al. ,Phys. Rev. D 66, 036001 (2002). 

[18] P. Mathews, Phys. Rev. D 63, 075007 (2001). 

[19] T. G. Rizzo Int. J. Mod. Phys. A18, 2797 (2003). 

[20] 1. Hinchliffe, N. Kersting, Y. L. Ma, Int. J. Mod. Phys. A19, 179 (2004). 

[21] M. R. Douglas and C. Hull, J. High Energy Physics 02, 008 (1998). 

[22] N. Seiberg and E. Witten, J. High Energy Physics 09, 032 (1999). 

[23] E. Witten, Nucl. Phys. B471, 135 (1996). 

[24] P. Horava and E. Witten, Nucl. Phys. B460, 506 (1996). 


14 


[25] A. Matusis, L. Susskind, and N. Toumbas, J. High Energy Physics 12, 002 (2000). 

[26] A. Bichl, J. Grimstrup, H. Grosse, L. Popp, M. Schweda, and R. Wulkenhaar, J. High Energy 
Physics 13, 0106 (2001). 

[27] G. P. Martin, Nucl. Phys. B652, 72 (2003). 

[28] Alvarez-Gaume L. and Vazquez-Mozo M.A., Nucl. Phys. B668, 293 (2003). 

[29] E. Brandt, G. P. Martin and E. Ruiz Ruiz, J. High Energy Physics 07, 068 (2003). 

[30] M. Buric, V. Radovanovic, and J. Trampetic, J. High Energy Physics 0703, 030 (2007). 

[31] Schupp P., You J., J. High Energy Physics 08, 107 (2008). 

[32] Horvat R., Trampetic J., J. High Energy Physics 1101, 112 (2011). 

[33] R. Horvat, A. Ilakovac, J. Trampetic and J. You, J. High Energy Physics 12, 081 (2011). 

[34] R. Horvat, A. Ilakovac, P. Schupp, J. Trampetic and J. You, [hep-th]]. J. High Energy Physics 
1204, 108 (2012). 

[35] R. Horvat, A. Ilakovac, J. Trampetic and J. You, J. High Energy Physics 1311, 071 (2013). 

[36] J. Trampetic and J. You, SIGMA 10, 054 (2014). 

[37] B. Jarco, L. Moller, S. Schraml, P. Schupp and J. Wess, Eur. Phys. J, 383 (2001). 

[38] T. Mehen and M. B. Wise, J. High Energy Physics 12, 008 (2000). 

[39] X. Galmet, B. Jurco, P. Schupp, J. Wess and M. Wohlgenannt, Eur. Phys. J. G 23, 363 (2002). 

[40] X. Galmet and M. Wohlgenannt Phys. Rev. D 68, 025016 (2003). 

[41] A. Prakash, A. Mitra and P. K. Das, Phys. Rev. D 82, 055020 (2010). 

[42] B. Melic, K. P. Kumericki, J. Trampetic, P. Schupp and M. Wohlgenannt Eur. Phys. J. G 42, 
483 (2005). 

[43] B. Melic et al , Eur. Phys. J. G 42, 499 (2005). 

[44] W. Behr, N. G. Deshpande, G. Duplancic, P. Schupp, J. Trampetic and J. Wess, Eur. Phys. 
J. G 29, 441 (2003). 

[45] P. Aschieri, B. Jurco, P. Schupp and J. Wess Nucl. Phys. B651, 45 (2003). 

[46] G. Duplancic, P. Schupp and J. Trampetic, Eur. Phys. J. G 32, 141 (2003). 

[47] M. Buric, D. Latas, V. Radovanovic and J. Trampetic, Phys. Rev. D 75, 097701 (2007). 

[48] M. M. Ettefaghi, M. Haghighat, Phys. Rev. D 75, 125002 (2007). 

[49] X. Galmet, Eur. Phys. J. G 41, 269 (2005). 

[50] X. Galmet and M. Selvaggi, Phys. Rev. D 74, 037901 (2006). 

[51] P. Schupp, J. Trampetic, J. Wess and G. Raffelt Eur. Phys. J. G 36, 405 (2004). 

[52] R. Horvat and J. Trampetic, Phys. Rev. D 79, 087701 (2009). 

[53] R. Horvat, D. Kekez and J. Trampetic, Phys. Rev. D 83, 065013 (2011). 

[54] Horvat R., Ilakovac A., Kekez D., Trampetic J., You J., J. Phys. G: Nucl. Part. Phys. 41, 
055007 (2014). 

[55] B. Melic et al , Phys. Rev. D 72, 054004 (2005). 

[56] G. Tamarit and J. Trampetic, Phys. Rev. D 79, 025020 ( 2009) 

[57] B. Melic, K. Passek-Kumericki and J. Trampetic, Phys. Rev. D 72, 057502 ( 2005) 


15 



[58] J. i. Kamoshita, Eur. Phys. J. C 52, 451 (2007). 

[59] M. Haghighat et al , Phys. Rev. D 82, 016007 (2010). 

[60] P. K. Das, N. G. Deshpande and G. Rajasekaran, Phys. Rev. D 77, 035010 (2008). 

[61] S. K. Garg, T. Shrecharan, P. K. Das, N. G. Deshpande and G. Rajasekaran, J. High Energy 
Physics 1107, 24 (2011). 

[62] Abhishodh, P, A, Mitra and Prasanta K. Das, Phys. Rev. D 82, 05020 (2010). 

[63] Prasanta K. Das and Abhishodh, P, Int. J. Mod. Phys. A27, 1250141 (2012). 

[64] Prasanta K. Das et al , Phys. Rev. D 83, 056002 (2011). 

[65] Prasanta K. Das and Abhishodh, P, Int. J. Mod. Phys. A28, 1350004 (2013). 

[66] W. Wang, E. Tian and Z. M. Sheng, Phys. Rev. D 84, 045012 (2011). 


16 



